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Abstract 

We point out that the proton-neutron energy contribution, for fow 
multipoles (in particular for the quadrupole component), effectively 
renormalizes the strength of the pairing interaction acting amongst 
identical nucleons filling up a single-j or a set of degenerate many- 
j shells. We carry out the calculation in lowest-order perturbation 
theory. We perform a study of this correction in various mass regions. 
These results may have implications for the use of pairing theory in 
medium-heavy nuclei and for the study of pairing energy corrections 
to the liquid drop model when studying nuclear masses. 
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1 Introduction 



The pairing force, which expresses in the most succint way the preference of 
nucleon pairs to become bound into .P = 0+ states in the atomic nucleus, 
has been widely used in many applications in the study of nuclear struc- 
ture properties ||, 0, ^. The special structure of the monopole pairing 
force has allowed to study the classification of nucleons that occupy a set of 
single-particle orbitals. The quasi-spin scheme 0, H, or, closely related the 
seniority quantum number v 0, leads to single-j and degenerate many-j shell 
exactly solvable models. They can be used as benchmarks to compare with 
approximation methods. Moreover, use of monopole pairing forces allows to 
determine an important energy contribution when evaluating total nuclear 
binding energies. 

Another important characteristic of the nucleon-nucleon effective inter- 
action acting inside atomic nuclei is expressed through the long-range com- 
ponents of this interaction [^]. The low multipoles and the quadrupole com- 
ponent, in particular, are essential in generating low- lying nuclear collective 
phenomena. They also contribute to the mean-field energy (binding energy 
in the nuclear ground state, deformation properties,...) of the atomic nucleus. 

These two components of the nucleon-nucleon effective force have formed 
a keystone to understand many facets of nuclear structure: from few valence 
nucleons near to closed-shell configurations as well as in those situations 
where many valence protons and neutrons are actively present outside closed 
shells. They are essential ingredients of any shell-model (or present large- 
scale shell-model) calculation, even though in most of them either model in- 
teractions are constructed explicitely or deduced from more realistic nucleon- 
nucleon potentials (see refs. 0, ^, ^, |Tl| to cite just some recent shell- 
model studies). 

It is our aim now, in the study of nuclear masses and two-neutron separa- 
tion energies S2n, to understand better the interplay of global energy contri- 
butions (liquid-drop energy as determined from simple models ITS] , [T^ or from 
more sophisticated microscopic-macroscopic methods |15], [16|, |18|, |T^) 
with local correlation effects. The latter contributions arise from specific 
nuclear structure effects taking pairing and low-multipole force components 
into account. Local correlation effects can come from various origins such 
as (i) the presence of closed-shell discontinuities, (ii) the appearance of local 
zones of nuclear deformation, and (iii) configuration mixing or shape mixing 
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that will show up in the ground state of the nucleus itself. A number of 
results have been published recently on this topic ^ . 

In the present paper, we study how monopole pairing (that forms an 
essential ingredient in all local energy correlations) can accomodate long- 
range forces and as such give rise to an effective pairing force that can later be 
used when (i) applying pairing theory in order to study lowest-order broken- 
pair excitations in medium-heavy nuclei, and (ii) see how, following results 
obtained recently by Fossion et al. |21| concerning the study of two-neutron 



separation energies S'2n, the new pairing corrections on top of the liquid-drop 
energy could reproduce local binding energy (and S2n) variations even better. 

In section 2, we succintly indicate the results of monopole pairing in a 
single-j shell, and we explicitely evaluate the proton- neutron quadrupole- 
quadrupole contribution to the ground-state energy. Thereby we observe 
that its effect leads to an effective pairing contribution. In section 3, we 
discuss applications in various mass regions in order to estimate the effect of 
this renormalized 'pairing-like'contribution. 

In the conclusion, we indicate that this extra effect may well be interesting 
when studying nuclear structure properties not too far from closed shells 
in which the interactions amongst identical nucleons dominate the proton- 
neutron interaction effects. 



2 Shell-model correlations 

As discussed in the introductory section, the interplay of the monopole pai- 
ring force and the proton-neutron low-multipole deformation-driving force 
components are essential ingredients of any shell model calculation. The 
aim of the present paper is to point out that, to lowest order, the proton- 
neutron quadrupole part can be incorporated as a renormalization of the 
pairing force strength. Of course, this means that applications will stick to 
regions near to closed shells where the proton-neutron energy contribution 
is not the dominant part. In open shells, with both valence protons and 
neutrons active, one has to treat both components (monopole pairing and 
low-multipole proton-neutron part) on equal footing. 

We first give a short reminder of the monopole pairing correlation energy 
part, considering the nucleons are filling a single-j shell (or a set of degenerate- 
j shells). Secondly, we study the proton-neutron energy correction (quadrupole 
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interaction) to be superposed to the monopole pairing part. 



2.1 Pairing energy corrections 

The Hamiltonian, describing the most simple case of a monopole pairing 
force between n identical valence nucleons interacting in a single-j shell, with 
a given strength G, is described as 

H = -G a]^a]_^a,_^,a,m'{-l?'^"'^'^' . (1) 

m,m'>0 

The ground state of such a system corresponds to a state where the nucleons 
are coupled in pairs of angular momentum J = 0. As a consequence the 
ground state has seniority t; = 0, 

|n,t; = 0) = (5;r/^|0), (2) 

where 

V" OT>0 

The binding energy of this system becomes 

^^pairing ^ '^(^2Vt - H + 2)n = G{n -N+ l)N, (4) 

where is the number of valence nucleon pairs and Q is the shell degeneracy, 
Q = j + 1/2. In a system with valence protons and neutrons, interacting 
through monopole pairing forces between alike nucleons, one has to consider 
the expression (|) for protons and for neutrons separately. As an illustration 
of this interaction (see figure |I|), we shown the spectrum of a monopole pairing 
force. From expression (Q), the two- neutron separation energy can be easily 
deduced, with result 

^™ = G{n, + 2-n,)= G{n, + 2 - 2N,), (5) 

where n^, is the number of valence neutrons, N^, the number of valence neu- 
trons pairs and Qi, = + 1/2 is the shell degeneracy for neutrons. 
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2.2 The quadrupole energy contribution 

As discussed before, we now evaluate the energy contribution in the ground 
state that is due to the proton-neutron quadrupole-quadrupole interaction. 
The quadrupole proton-neutron Hamiltonian can be written as, 

H = kQtt ■ Qu, (6) 

where Q is the quadrupole operator for protons and neutrons, respectively, 
and ■ stands for the scalar product. The main characteristic of the quadrupole 
operator is that it induces the breaking of pairs, promoting pairs coupled to 
J'^ = O"*" into pairs coupled to J'^ = 2~^, thereby changing the seniority 



quantum number v and causing a core-polarization effect [^. To a first 
approximation, the ground state of the system will change from a condensate 
of proton and neutron pairs coupled to J'^ = 0"^ into a superposition of the 
state of expression together with a new one where one-proton and one- 
neutron pairs coupled to J'^ = 2+ are induced, 

® N,) = \Sf;' ® Sfj; J = 0) + ^\Sf^--'Dj^ ® S^r^D,,- J = 0),(7) 

where the operator 

= E i-iy^'^ijmj - m\20)a]^a]_^, (8) 



will create a pair of nucleons coupled to J = 2. 

Taking into account that the quadrupole energy contribution is small with 
respect to the monopole pairing interaction energy, the most straightforward 
way in order to fully determine the state vector in expression (|^ is to use 
perturbation theory. The mixing coefficient results as 

where we use the shorthand notation 

a = {Sf^- ® Sfj'; J = 0\Q^ ■ Q,\Sfj''D,^ ® Sfr'D,,- J = 0), (10) 
and define the energy difference 

AE„ = -E2+-^2+, (11) 
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in which E2+ (p = vr, u) denotes the excitation energies of states with one 
J = 2 proton or neutron pair, respectively. 

The energy correction due to the quadrupole force is then given by eval- 
uating the matrix element of the Hamiltonian (^) using the state vector (|^). 
For a system where the forces are monopole pairing and quadrupole proton- 
neutron interactions only, this results into the binding energy expression 

BE = 5EP^™s + 2—— -, (12) 

{—AEa) 

where _BiJP'^'''"^s corresponds to the result given in expression 

Due to the schematic structure of the states appearing in expression ([10|), 
it is now possible to obtain an explicit expression for this matrix element [^]. 
In a first step, we decouple the proton and neutron parts with the result 

X {{S,,f-,J = 0\\Q4{Sj,f^-'D,,;J = 2). (13) 
Using the appropriate reduction formulae one arrives to the result 

a = -L^N^{n^-N^)^N,in,-N,) . ^^-^^^-^ (14) 
where, 

m = (SjAqadj.) = (15) 

and 

(Q.) = {s.AQAD,.) = -^=^=={UQAj.)- (16) 

The expressions (|15D and (0) can be even more simplified if one uses har- 
monic oscillator wave functions and uses the expression for the quadrupole 
operator, 

Q = T.riUO^,^^), (17) 
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resulting into 



3 

-{N,o + h , (18) 
2 ^(2j,-l)j,(j, + l)(2j, + 3) 

Here p = (tt, z/), l2(^, V^) denotes the spherical harmonic with L = 2 and 
describes the number of quanta of the shell. 

One finally obtains a closed expression for the binding energy, 

If we are interested in the study of binding energies within a set of isotopes 
(thus fi^r and are fixed numbers), and by defining the coefficient C as 
follows, 

one obtains a "correlated" binding energy expression in the ground state 

^^corrclatcd _ ^^pairing ^ ^^^^ - A^^ + 1)N^ - CN^. (21) 

We can incorporate the proton-neutron quadrupole binding energy contribu- 
tion as a renormalisation of the pairing strength in regions near closed shells 
where pairing dominates the proton-neutron quadrupole interaction. In these 
regions and for medium-heavy nuclei, which can be characterized by single- 
particle orbitals with large degeneracies fli, (see also discussion in sect. 3 and 
tables 1 and 3), the energy contribution for the third term in eq. (^Tj) is much 
smaller than the energy contribution for the second term at the beginning 
of the shell (small N^) and up towards midshell {N^, ^ The total 

"correlated" binding energy then becomes to a good approximation, 

^^corrclatcd ~ (G- + C) {Q^ - + l)N^. (22) 

So, one obtains a form, identical to the original monopole pairing energy 
expression, albeit with a new coupling strength. This equation can be applied 
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to the region where detailed values for binding energies in the Z ~ 40 and 
Z ~ 50 nuclei (see table 3) are known. 

In the next section we shall evaluate, in some detail, the range of values for 
the coefficients G and C in medium-heavy nuclei. These results will show how 
good our idea is in reality, that proton-neutron quadrupole-quadrupole forces 
can be used to define an effective pairing force between identical nucleons. 

An interesting result is that in evaluating the two-neutron separation 
energy, the quadratic part drops out and one obtains a strictly linear behavior 
in Ny with the result 

^quadrupole^^^^ = (G + (fi, + 2 - - C. (23) 

This is essentially the same result as was obtained using a pure monopole 
pairing force (see expression (^) except for the small correction factor C 
(which we still have to prove). 

The former discussion which uses first order perturbation theory to de- 
termine the wave function (eq. (|^)) can be repeated for pure hexadecupole or 
higher multipole interaction Hamiltonians (using eq. (|^), replacing Qp by the 
appropriate multipole operator). Each higher multipole separately will give 
a smaller energy contribution to the binding energy, so that only a limited 
number of multipoles will be important. Using higher order perturbation 
theory to determine the modified wave functions, extra energy contributions 
become are possible that come from different multipoles acting together. 
These higher order contributions result in a different dependence than 
the contributions from separately treated multipoles. These higher order ef- 
fects will be of minor importance. In the present paper, we do not aim at 
carrying out a detailed shell model study - in that case the multipole Hamil- 
tonian should rather be diagonalised - but we find it quite surprising that the 
lowest-order effect induced by proton-neutron residual forces comes about as 
a 'renormalization' of the original monopole pairing force. We find this an 
interesting observation. 
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3 The pairing and quadrupole strength: some 
specific apphcations and how well works 
the above approximation 

In the present section, we make a detailed study of the correction factor 
C , mainly concentrating on medium-heavy nuclei. We compare this effec- 
tive pairing force with the monopole pairing strength that is derived from 
standard parametrizations, in this mass region. 

• The monopole pairing strength G. 

An average value of G for medium-mass and heavy nuclei is 25/A MeV 



23| . However, this value can significantly change in different mass 
regions, ranging from G = 19/A MeV, for the regions Z ^ 40 — 50, to 
^ = 30/^ MeV for Pb nuclei [H. 

As an example, for 52Te isotopes {A ^ 106 — 152) the pairing strength 
is G ^ 0.15 MeV, and for 42 Mo isotopes {A ^ 84-128) G ^ 0.18 MeV. 

The effective pairing strength G, obtained from the proton-neutron 
quadrupole force. 

In order to estimate the value of G, using the method discussed in 
section 2, we have to reduce the set of more realistic single-particle 
orbitals that are not degenerate and are typical for a given mass region, 
into a large, degenerate single-j shell. In table |I|, we indicate those 
values of j as well as the degeneracies for different major shells or 
subshells in the regions Z ~ 40 and Z ~ 50. As example of isotopes 
with Z ~ 40 we consider the 48Cd (N^ = 1), 42 Mo (N^, = 1), 44 Ru 
(N^ = 2), and 46Pd (N^ = 2) nuclei . For the region with Z ~ 50, we 
consider the 52 Te (N^ = 1), 54Xe (N^ = 2), seBa (N^ = 3), and sgCe 
(Ntt = 4) nuclei. For both shells we take as number of quanta = 4, 
which covers a shell from 40 to 70 nucleons and is in agreement with 
table m. 

Another ingredient necessary in order to determine G (see expression 
(|20|) ) is the value of AEa = —E2+ — E2+ which can be estimated from 
experimental energy systematics of the first 2"*" excitation energy [^]. 
This results for the Z ~ 50 region into a maximum value of AE^ — —2 
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MeV and for the Z ~ 40 region, a slightly higher value /S.Ea — —2.5 



MeV is obtained 22 



The final element in order to calculate C is the strength of the proton- 
neutron interaction k. There exist realistic values of the quadrupole 



strength in the context of the IBM [26, 27], kq, that can be easily related 



with the value of k appearing in expressions (B|) and (EDF) through the 



following relation |22]: 



As can be verified from table and noted in figure |^, the reduced matrix 
elements |(jp||<5p||jp)| increase linearly with jp, but this variation is 
almost completely compensated through the presence in the numerator 
of the factors Vtp. This then results in values for k that are almost 
independent of the value of jp iytp). On the other hand, the values 
of {Qp) and {Qp)"^ (the third factor in expression (pOD ) do not change 
dramatically with jp in the range Jp = | — This is an important 
outcome which implies that the value of C is not strongly dependent 
on the particular j value of the proton and/or neutron orbitals 
(degeneracies) that the nucleons are occupying (state independence). 
This results into robust values of C. 

Having discussed the various elements that are necessary in order to 
derive a schematic, albeit rather realistic, estimate of C, we present, 
in table ^, the values for different nuclei in the regions Z ~ 40 and 
Z ~ 50. We point out that two alternatives shell closures have been 
used. 

Inspecting the results, as given in table |^, it becomes clear that the 
effective pairing correction, stemming from the proton-neutron interac- 
tion, can be at maximal of the order of 10-15% of the regular monopole 
pairing strength. The precise values depend somewhat on the degenera- 
cies of the proton and neutron shell-model spaces that are used in order 
to describe the various series of isotopes. 

As a conclusion, on notices that the quadrupole constant C is smaller 
than the monopole pairing constant G by a factor varying in between 5 to 
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25. Therefore, quadrupole and low-mult ipole force components give rise to 
contributions to the binding energy that exhibit the same structure (except 
for some smaller corrections) as those resulting from monopole pairing forces 
solely. The particular number dependence of the binding energy (see ex- 
pression (p2D) occurs through the common factor {Q — N + 1)N with Q the 
degenarcy of the shell that is filling up with identical nucleons. 



4 Conclusion 

In the present paper, we have shown that pure monopole pairing can acco- 
modate long-range forces (we have studied in particular the case of proton- 
neutron quadrupole forces but the extension to other low multipoles is now 
straightforward) using perturbation theory, and as such give rise to an ef- 
fective pairing force. This latter effect renormalizes the monopole pairing 
energy with an amount that can vary between 5 to 25%. 

This result came partly as a surprise and did come up when we were 
studying the variation of nuclear binding energies and, more in particular, 
two-neutron separation energies 5*2™ over a large region of nuclei (rare-earth 



mass region, nuclei in the neutron-deficient Pb region |20, Pl[). In the above 



papers, we have studied the local correlation energy within the framework of 



the Interacting Boson Model (IBM) p6|, |28| and outlined a prescription 
in order to derive nuclear masses within a single framework. 

Recently, we have started the study of nuclear binding energies and two- 
neutron separation energies in medium-heavy nuclei (mass A=100-130 
region) using the same concepts but also trying to evaluate the local corre- 
lation energy taking into account the shell-model structure and the residual 
interactions (pairing, quadrupole proton- neutron interactions) |2^. In order 
to study these local energy corrections that appear on top of the global energy 
(liquid drop energy, macroscopic-microscopic energy studies - see refs. given 
in the introduction), we have studied in the present paper the modifications 
that proton-neutron forces induce on the strict monopole pairing energy. We 
have shown that the dependence on nucleon number (e.g. neutron number 
when studying series of isotopes) for the effective pairing contribution is iden- 
tical to the nucleon number dependence of the monopole pairing force. In 
going through a series of isotopes, changing the number of neutron pairs, N^,, 
the effective pairing force will likewise exhibit a neutron number dependence 
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through the presence of the factor C (see expression pO|). 

Therefore, we aim at (i) applying pairing theory to study lowest-order 
broken-pair excitations in medium-heavy nuclei, and (ii) see how, following 
up on results obtained recently by Fossion et al. [|21[ concerning the study 
of two- neutron separation energies S'2„, the effective pairing corrections on 
top of the liquid-drop energy could reproduce local binding energy (and 82-0) 
variations even better. 
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Table Captions 



Table 1: Degeneracies and associated single-j shell configurations for regions 
Z ~ 40 and Z ~ 50. 



Z~ 40 


Shell 


Degeneracy 


j 


Z = 40 - 50 
iV = 50 - 64 
= 50 - 82 


n, = 7 
= 16 


Jn = 9/2 
Ju = 13/2 
J. = 31/2 


Z ~ 50 


Shell 


Degeneracy 


J 


Z = 50 - 82 

Z = 50 - 64 
= 50 - 82 


= 16 
= 16 


jn = 31/2 

J. = 13/2 
j. = 31/2 
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Table 2: Values for (Qp) , (Qp)^ and {jp \\Qp\\ jp) (see expressions |T^, |T^ and 
p!8|) for jp going from 3/2 to 31/2. 



Jp 




{QpY 


{JpWQpWJp) 


3/2 


-3.10 


9.62 


-3.10 


5/2 


-3.31 


11.00 


-4.06 


7/2 


-3.38 


11.46 


-4.78 


9/2 


-3.41 


11.67 


-5.40 


11/2 


-3.43 


11.78 


-5.94 


13/2 


-3.44 


11.85 


-6.44 


15/2 


-3.44 


11.89 


-6.89 


17/2 


-3.45 


11.92 


-7.32 


19/2 


-3.45 


11.94 


-7.72 


21/2 


-3.45 


11.96 


-8.11 


23/2 


-3.46 


11.97 


-8.47 


25/2 


-3.46 


11.98 


-8.82 


27/2 


-3.46 


11.99 


-9.48 


31/2 


-3.46 


12.00 


-9.79 
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Table 3: Values of C ( in MeV) for different nuclei and different mass regions. 



Z ~ 40 


Shells 


Mo 


Cd 


Ru 


Pd 


Qtt = 5, fli, = 16 

= 5. n,. = 7 


0.0062 

0.0059 


0.0062 

0.0059 


0.0093 

0.0079 


0.0093 

0.0079 


Z ~ 50 


Shells 


Te 


Xe 


Ba 


Ce 


= 16, = 16 

= 7 , fi,^ = 16 


0.0091 
0.0036 


0.017 
0.006 


0.0240 
0.0073 


0.029 
0.0073 
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Figure Captions 



Figure 1: Spectrum for a pure pairing force within the 1/111/2 orbital, i.e. for 
the (1/^11/2)" spectrum, with seniority v = 0, v = 2, v = 4, and = 6 as a 
function of the particle number n. The pairing strength G has been chosen 
as G = 0.25 MeV M. 
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Figure 2: Plot of the values for (Qp) and {jp \\Qp\\jp) for jp going from 3/2 
to 31/2. As one can see, {Qp) is nearly a constant for changing jp. 
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